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ABSTRACT 

Values of the Hubble constant reported to date which are based on measurement of 
the Sunyaev-Zeldovich (SZ) effect in clusters of galaxies are systematically lower than 
those derived by other methods (e.g., Cepheid variable stars, or the Tully-Fisher rela- 
tion). We investigate the possibility that systematic errors may be introduced into the 
analysis by the generally adopted assumptions that observed clusters are in hydrostatic 
equilibrium, are spherically symmetric, and are isothermal. We construct self-consistent 
theoretical models of merging clusters of galaxies using hydrodynamical/N-body simu- 
lations. We then compute the magnitude of H D derived from the SZ effect at different 
times and at different projection angles both from first principles, and by applying each 

iCurrent Address: Goddard Space Flight Center, Code 930, Greenbelt, MD 20771 



1 



of the standard assumptions used in the interpretation of observations. Our results indi- 
cate that the assumption of isothermality in the evolving clusters can result in H Q being 
underestimated by 10-30% depending on both epoch and projection angle. Moreover, 
use of the projected, emission-weighted temperature profile under the assumption of 
spherical symmetry does not significantly improve the situation except in the case of 
more extreme mergers (i.e., those involving relatively gas-rich subclusters). Although 
less significant, we find that asphericity in the gas density can also result in a 15% error 
in H D . If the cluster is prolate (as is generally the case for on- axis, or nearly on-axis 
mergers), and viewed along its major axis, H will be systematically underestimated. 
More extreme offaxis mergers may result in oblate merger remnants which when viewed 
nearly face-on may result in an overestimation of H Q . A similar effect is noted when 
viewing a prolate distribution along a line-of-sight which is nearly perpendicular to its 
major axis. In both cases the potential overestimation occurs only when the remnant is 
viewed within 15-30° of face-on. Bulk gas motions and the kinematic SZ effect do not 
appear to be significant except for a brief period during the very early stages of a merger. 
Our study shows that the most meaningful SZ measurement will be accompanied by a 
high resolution temperature data and a detailed dynamical modeling of the observed 
system. In lieu of this, a large sample selected to avoid dynamically evolving systems is 
preferred. 

Subject headings: Galaxies: clusters: general - galaxies: intergalactic medium - 
hydrodynamics - methods: numerical 
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1. INTRODUCTION 

The Hubble constant (H ) and determination of 
its value are among the most controversial topics 
within the astronomical community. Recently, how- 
ever, there seems to be a general convergence toward a 
value for H at the high-end of the currently accepted 
range, (i.e., > 70 km s _1 Mpc -1 ). Such a value is sup- 
ported by a variety of methods including the Tully- 
Fisher relation (85±10; Pierce & Tully 1988), galaxy 
surface brightness fluctuations (82±7; Tonry 1991), 
Type II supernovae and the expanding photosphere 
method (EPM) (73±7; Schmidt et al. 1994), as well 
as both space-based (80±17; Freedman et al. 1994) 
and ground-based (87±7; Pierce et al. 1994) mea- 
surements of Cepheid variable stars in the Virgo clus- 
ter. Mould et al. (1995) provides an excellent review 
of the recent analysis applied to the Virgo cluster in 
which they note that all of the above numbers are 
consistent with 80±17 km s _1 Mpc -1 . 

In contrast to these values, two methods have con- 
sistently supported a lower (< 60 km s -1 Mpc -1 ) 
value for H . Schaefer (1996) reviews the observations 
of 10 Type la supernovae (SNIa) and reports a mean 
value of 55±3. However, Hamuy et al. (1995) finds 
that values of H based on SNIa will be systemat- 
ically underestimated by 15% unless a correction is 
applied using the luminosity-decline rate relation of 
Phillips (1993) which leads to a best value near 65 
km s -1 Mpc -1 . This result is consistent with the 
SNIa-based value derived by Riess, Press, & Kirshner 
(1994), 67±7kms -1 Mpc -1 . Mould et al. (1995), us- 
ing a recalibration of SNIa absolute magnitudes based 
on six SNIa in the Virgo cluster, find a value of 71±7 
km s -1 Mpc -1 . 

The second method which has consistently sup- 
ported lower values of H , and the one of interest here, 
relies on the measurement of the Sunyaev-Zeldovich 
(SZ) effect (Sunyaev & Zeldovich 1972). For a com- 
plete and recent review of the SZ effect, the reader 
is referred to Rephaeli (1995). Briefly, the SZ effect 
is a distortion of the cosmic microwave background 
(CMB) spectrum caused by inverse Compton scatter- 
ing of low energy CMB photons by high energy elec- 
trons in the dense cores of massive clusters of galaxies. 
Knowing the density and temperature of the scatter- 
ing plasma, as determined from X-ray observations of 
the intracluster medium (ICM) , and the CMB bright- 
ness temperature in the cluster core, as measured by 
millimeter observations, it is possible to determine 



H (see sQ). Owing to the difficulty of the mil- 
limeter observations and, until recently, the limited 
quality of X-ray temperature data, there have only 
been a few determinations of H using this method. 
Still, they have tended to support a lower value for 
H (e.g., 45±17, Birkinshaw et al. 1991; 48±28, Jones 
et al. 1993; 55±17 Birkenshaw & Hughes 1994). The 
one notable exception is a recent measurement using 
the Coma cluster which finds H Q = 71±25 km s -1 
Mpc -1 (Herbig et al. 1995). This value is based 
on a Ginga temperature for Coma of 9.1 keV. New 
ASCA observations resulted in a central temperature 
for Coma of ~8 keV (Honda et al. 1996) which is 
more consistent with previous temperature estimates 
such as 7.5±0.2 keV from Tenma, 8.5±0.3 keV from 
EXOSAT (Hughes 1989), and 7.7 keV from Spacelab 
2 (Watt et al. 1992). Using the lower temperature 
estimate would reduce the the Coma-derived value of 
H to a value nearer 50 km s -1 Mpc -1 . The review ar- 
ticle by Rephaeli (1995) lists two other measurements 
(A2256 and A2163) of H using the SZ effect which 
are greater than 70 km s -1 Mpc -1 . However, both 
of these cluster show significant substructure, neither 
of these H measurements have appeared, as of yet, 
in the rcfcrccd journals. Similarly, Kobayashi, Sasaki 
& Suto (1996) argue that SZ based measurements are 
consistent with those of Cepheids. This consistency is 
owed to the large uncertainty in the SZ derived values 
of H . 

What are the sources of this apparent systematic 
discrepancy? Although the X-ray data have improved 
dramatically in the last decade, it is still difficult to 
determine the internal structure of clusters from the 
X-ray imaging because it only supplies projected tem- 
perature and surface brightness information, while 
supplying no information regarding the internal gas 
dynamics. Consequently, observers are often forced 
to assume that clusters are spherical, hydrostatic, and 
isothermal systems. Such assumptions are contrary to 
recent observational evidence that shows many clus- 
ters (30-70%, Jones & Forman 1991; Bird 1993; Davis 
1994; Mohr et al. 1995) are still dynamically evolv- 
ing. In a few instances, detailed temperature maps of 
clusters have been produced only to reveal extremely 
complex temperature morphologies (e.g., 2256, Briel 
& Henry 1994; Roettiger, Burns & Pinkney 1995; 
A754, Henriksen & Markevitch 1996) We attempt to 
determine the source of the discrepancy using self- 
consistent numerical hydrodynamical/N-body simu- 
lations of evolving clusters of galaxies. Specifically, 
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we address systematic errors associated with the in- 
herent limitations in the X-ray data. Our basic ap- 
proach is to "observe" massive simulated clusters at 
several epochs during mergers with a variety of sub- 
clusters. We focus on mergers, because it is believed 
that massive clusters grow through the accretion of 
smaller structures, and mergers are the only mech- 
anism likely to affect a cluster's global properties. 
Using the same techniques commonly employed by 
observers, we determine the temperature and den- 
sity profiles of the simulated clusters as a function of 
viewing angle and epoch. We then compare the "ob- 
served" value of H with the "true" value calculated 
from the line-of-sight (LOS) temperature and density 
profiles, i.e., from first principles. By comparison of 
the observed and true values of H , we can evaluate 
the systematic errors. 

It should be emphasized that we are analyzing the 
remnants of strong mergers which are therefore dis- 
rupted systems, and in some sense, they may repre- 
sent a worst case scenario. On the other hand, it 
should be kept in mind that selection effects exist 
such that a sample of clusters having SZ measure- 
ments can potentially be biased toward evolving sys- 
tems such as these. We find that such systems are 
hotter and, consequently, more X-ray luminous then 
their relaxed counterparts. Both effects will bias their 
selection. In addition, we find merger remnants to be 
prolate systems. Consequently, they will have an en- 
hanced X-ray surface brightness when viewed end-on 
which may further bias their selection. It should be 
noted that all but two of the clusters listed in the 
sample presented by Rephaeli (1995) contain signifi- 
cant substructure, and those two (A478 and A2142) 
are considered to be strong cooling flows (Edge et 
al. 1992), a matter that further complicates the SZ 
analysis. 

In §^|, we describe the Sunyaev-Zeldovich effect 
and how it is used to determine H . Section |3] de- 
scribes our numerical methodology and initial con- 
ditions for computing cluster mergers, as well as, our 
choice of parameter space, our data analysis, and tests 
of our method. Section [h] contains a physical descrip- 
tion of the cluster mergers. Section ^ contains our 
major results and a comparison with results from a 
previous study by Inagaki et al. (1995). In §^], we 
discuss the origin of the observed discrepancy, limi- 
tations of this study, and implications for observers. 
And, finally, we summarize our conclusions in §(H 



2. THE SUNYAEV-ZELDOVICH EFFECT 
AND THE HUBBLE CONSTANT 

The SZ effect is a spectral distortion of the cos- 
mic microwave background (CMB) resulting from in- 
verse Compton scattering of low energy CMB pho- 
tons by high energy electrons in the ICM (Sunyaev & 
Zeldovich 1972). The change in the CMB brightness 
temperature observed in the Rayleigh- Jeans region of 
the spectrum is: 

= -2^4 / n M (1) 
1 CMB m e c z J 

where n e is the electron density, T e is the electron 
temperature, ot is the cross section for Thomson 
scattering, and the line integral, dl , is performed 
along the LOS through the cluster. Since n e oc h 1 / 2 
and dl oc h^ 1 (T e is independent of h), we can solve 
Eq. [| for h in terms of the observable quantities AT 
(millimeter observations), n e {r) and T e (r) (X-ray ob- 
servations) , 

A more detailed derivation of the dependence of 
H D on these quantities can be found in Birkinshaw 
et al. (1991) and Rephaeli (1995). Rephaeli (1995) 
also includes the relativistically correct form. 

Observationally, the situation is complicated by 
several factors. Typical values for the temperature 
decrement are of order a few hundred fiK thus re- 
quiring extreme sensitivity of the millimeter measure- 
ments which in turn limits this type of observation 
to the hottest, densest, and presumably most mas- 
sive clusters. Often, these are distant clusters which 
are affected by low spatial resolution in both the mil- 
limeter and X-ray observations. Although low spa- 
tial resolution and the use of a reference beam in 
the millimeter observations introduces a significant 
model dependence to the analysis, the central av- 
eraging may actually help smooth small scale tem- 
perature fluctuations and thereby help mitigate some 
of the associated uncertainties (Herbig 1996). The 
typical resolution of the millimeter observations are 
>V (e.g., 1.78' , Birkinshaw & Hughes 1994; 7.3' , 
Herbig et al. 1995). It is therefore necessary to make 
assumptions about the cluster gas distribution in or- 
der to deduce the central temperature decrement from 
the observed beam averaged value (Birkinshaw et 
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al. 1984). Only recently have interferometric tech- 
niques allowed the generation of SZ maps with sub- 
arcmin resolution (e.g., Carlstrom et al. 1996). In 
addition to these sources of error, observers must con- 
tend with contamination by cluster radio sources. 

In this study, we do not address the uncertainties 
associated with the measurement of AT/T. There- 
fore, from Eq. |^, the fractional error in H based 
solely on limitations of the X-ray data is simply the 
square of the ratio of the observationally derived pres- 
sure integral to the true LOS pressure integral or, 



h , 



■bs 



[Jn e T e dl] o 
[J n e T e dl] tr 



(3) 



where the RHS numerator is based on observationally 
derived values for n e , T ei and dl and the denominator 
is based on the true 3-dimensional simulation data. 
Both integrals are calculated at the resolution of the 
simulation, 50 kpc (See §||]). 

3. METHOD 

3.1. Dynamical Calculations of Cluster Merg- 



3.1.1. Numerical Method 

To compute the dynamical evolution of merging 
clusters of galaxies, we use a hybrid hydrodynamical/N- 
body code in which the hydrodynamical component is 
CMHOG written by one of us, J. M. Stone. CMHOG 
solves the fluid equations using an implementation 
of the piecewise-parabolic method (PPM, Colella & 
Woodward 1984) in its Lagrangian remap formula- 
tion for gas dynamics. The primary advantages of 
PPM are lower numerical diffusion through the use 
of third order parabolic spatial interpolations, and its 
ability to accurately resolve shocks in as few as 1- 
2 zones through the use of a Ricmann solver. The 
result is an increase in the small structure captured 
at a given resolution over codes employing artificial 
viscosity to capture shocks. This code has been ap- 
plied to a variety of astrophysical problems and has 
been extensively tested using the problems described 
in Woodward and Colella (1984). 

The collisionless dark matter is evolved using an 
N-body code based on a standard particle-mesh algo- 
rithm (PM). The particles are evolved on the same 
grid as the gas using the same time step. The time 
step is determined by applying the Courant condition 



simultaneously to both the dark matter and the hy- 
drodynamics. The only interaction between the colli- 
sionless particles and the gas is via the gravitational 
potential. Therefore, the two codes are run largely 
independent of each other. The hydrocode supplies 
the gas density distribution to the N-body code as 
part of the source function to Poisson's equation (the 
self-gravity of both the gas and dark matter is in- 
cluded in these simulations). The N-body code, in 
turn, supplies the gravitational potential to the hy- 
drocode where it becomes a source function to the 
energy equation. Since we are modeling an isolated 
region, the boundary conditions for Poisson's equa- 
tion are determined by a multipole expansion of the 
mass distribution contained within the grid. Particles 
that leave the grid are lost to the simulation. Typi- 
cally, less than a few percent of the particles leave the 
grid. 

The computational grid is similar to that described 
in Roettiger, Loken, & Burns (1997). The simulation 
is fully three-dimensional. The grid is a fixed and 
rectangular (100 x 100 x 220 zones). The merger axis 
coincides with the grid's major axis along which res- 
olution is uniform and scales to 50 kpc or 6 zones 
per primary cluster core radius. Resolution along 
the grid's minor axis is uniform within the central 
40 zones and ratioed in the 30 zones on either side. 
That is, the resolution is a uniform 50 kpc extend- 
ing to 1 Mpc (20 zones) on either side of the merger 
axis. Beyond 1 Mpc, the zone dimensions increase by 
3% from one zone to the next out to the edge of the 
grid. We use outflow boundary conditions for the hy- 
drodynamical evolution. We do not include radiative 
cooling in these simulations. 

3.1.2. Initial Conditions 

The initial conditions are similar to those in sev- 
eral previous studies (Roettiger, Loken & Burns 1993; 
Roettiger, Burns & Pinkney 1995; Roettiger et 
al. 1996, 1997). We begin with two clusters that 
are initially isothermal (r < 4r c ) and in hydrostatic 
equilibrium. They are simply placed on the com- 
putational grid separated by ~ 5 Mpc and allowed 
to merge under the influence of their mutual grav- 
ity. The clusters are given an initial relative velocity 
of 300 km s _1 . While not significantly affecting the 
final impact velocity, the small initial velocity does 
greatly reduce the time required to merge and thus 
saves considerable computational resources. The ide- 
alized initial conditions have several benefits. They 
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allow a more efficient use of the computational vol- 
ume, and consequently, higher resolution. They pro- 
vide a well-defined baseline with which to compare 
the subsequent evolution of the post-merger clusters. 
And, they allow for a systematic study of merger pa- 
rameter space. 

Whereas the basic setup is identical to previous 
simulations by Roettiger et al. , the details of the 
intial cluster structure have changed somewhat. The 
total mass distribution in these simulations is well- 
characterized by a power-law slope of approximately 
-2.5 for r > r c . The gas distribution is typically some- 
what flatter, ranging from -2.0 to -2.5 depending on 
the intial (3 parameter, see Table 1. The tempera- 
ture profiles are also somewhat different from previ- 
ous simulations. These clusters are strictly isother- 
mal only within 4r c . At larger radii, the tempera- 
ture decreases gradually. Figure [l] shows an example 
of the initial density and temperature profiles for a 
P~l cluster. Here, we define two values of /?, des- 
ignated Pfu and pspec- The value of /3fu is deter- 
mined by fitting an isothermal /3 model to the X- 
ray surface brightness distribution (see eq. 5), while 
Pspec = ^m p a v 2 /kT where /i is the mean molecular 
weight, Trip is the proton mass, k is Boltzmann's con- 
stant, a v is the galaxy velocity dispersion, and T is 
the ICM temperature. Initially, Pfu= Pspec, at least 
within the central 3-4r c . Unless otherwise stated, all 
references to p refer to ft fit- Finally, since these sim- 
ulations are non-cosmological, the value of H does 
not directly influence the scaling of the cluster di- 
mensions. However, the arbitrary length scaling that 
we have chosen most closely resembles clusters in a 
universe m which H =50 km s" 1 Mpc" 1 . 

3.1.3. Parameter Space 

In order to constrain our merger parameter space, 
we need to ask several questions. First, given a mas- 
sive cluster of the type necessary for the SZ analysis, 
what mass of cluster is it most likely to have recently 
interacted with it? We gain some insight into this 
matter by looking at the mass function for clusters. 
Bahcall & Cen (1993) analytically represent the clus- 
ter mass function as, 

n(> M) = 4xl0- 5 (M/M*)- 1 exp(~M/M*)h 3 Mpc- 3 

(4) 

where n(> M) is the number density of clusters with 
mass greater than M, M* = 1.8 ± 0.3 x lO 14 /^ 1 M Q , 
and h = H o /100. Normalizing to a region of space 



large enough to have one cluster of mass comparable 
to our primary (~10 15 M ), we find that -90% of 
the other clusters in this same region are <25% of 
its mass. Thus, we focus on mergers between clusters 
with large mass ratios. 

A second consideration is the relative gas content 
of the two clusters, particularly their central gas den- 
sities which, when combined with the impact velocity 
determines the ram pressure experienced by the re- 
spective clusters. Previous simulations have shown 
the importance of the self-interaction of the gas com- 
ponents in determining the post-merger evolution of 
the ICM. The relative gas densities will determine the 
degree to which the subcluster gas is stripped, the 
strength of shocks generated during the merger, and 
the accompanying heating, as well as the size and du- 
ration of bulk flows and the rate at which the remnant 
returns to equilibrium. X-ray studies of clusters of 
galaxies show the typical range of central gas densities 
to be 10~ 4 < n < 10~ 2 cm" 3 (Sarazin 1986), while 
the global gas content of clusters, by mass, seems to 
range from a few percent to >30% with typical values 
in the 10-20% range (White & Fabian 1995). 

Since we are modeling a wide range in total clus- 
ter mass (up to a factor of 8), it would be helpful to 
constrain the models if there were a correlation be- 
tween total mass and the ICM properties. Assuming 
hydrostatic equilibrium is the norm rather than the 
exception, one would expect a reasonable correlation 
between total mass and ICM temperature. However, 
a similar correlation does not appear to exist between 
total mass and gas density or baryon fraction (White 
& Fabian 1995). In the absence of some extreme local 
physics (e.g., AGN, star formation), one would expect 
the fraction of mass in baryons to be relatively con- 
stant from one cluster to the next. The fact that we 
observe a wide range of values is, in and of itself, a 
very interesting problem with considerable cosmolog- 
ical implications (e.g., Lubin et al. 1996). Still, it 
does little to help constrain our models. The situa- 
tion is no better with the central gas density. Here, 
studies are limited by the resolution of the X-ray ob- 
servations, and it may be more reasonable to expect 
local physics, particularly radiative cooling, to play a 
significant role. It has been noted that many very low 
mass systems tend to have more extended gas distri- 
butions than their high mass counter parts (i.e., low 
P) (Doe et al. 1995; Mulchaey et al. 1996). In which 
case, if the gas contents are comparable and the dark 
matter distributions similar, one might then expect 
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at least on average, that lower mass systems will have 
lower central densities. 

Table 1 contains the parameters of the clusters 
used in this analysis. In light of the above consid- 
erations, we have run six head-on merger simulations 
in which we vary the relative cluster masses and gas 
properties. In each case, the primary cluster is scaled 
to 1.5 x 10 15 M Q . We have evolved both 4:1 and 
8:1 total mass ratio mergers. Within each mass ratio, 
we vary both the gas content and central gas densi- 
ties within the observed range of parameters. Central 
gas density varies from ~8x 10~ 4 to nearly 2 xlO -2 
cm~ 3 . The fraction of mass in baryons within the cen- 
tral 1 Mpc ranges from 0.08 to 0.22 with most clus- 
ters ^0.12. Merger 7 is an offaxis merge r in which 
the impact parameter is ~160 kpc (see §5^5.). Al- 
though this study in no way represents a complete 
sampling of parameters space, wc do believe it to be 
representative of that space most likely to affect the 
SZ analysis. Mergers with lower mass subclusters or 
those with lower central gas densities relative to the 
primary are less likely to significantly influence the 
internal structure of the primary cluster. 

3.2. Synthetic Data Analysis of Merger Sim- 
ulations 

In this section we describe the manner in which we 
attempt to obtain the X-ray observables, T e and n e , 
from the simulated data base. It is our goal to mimic 
the observational procedures as closely as possible in 
order to examine the effects of limited information in 
the observational data. We do not address problems 
associated with observational concerns such as photon 
statistics or poor spectral resolution. We therefore as- 
sume the clusters are well-observed within the region 
being analyzed. 

3.2.1. Observed Temperature 

Observationally, we are interested in the projected, 
emission-weighted ICM temperature. After selecting 
a particular epoch within a given merger, we calcu- 
late the X-ray volume emissivity (ey) within a chosen 
bandpass (e.g., ASCA, 0.5-10 keV). We use the Mewe- 
Kaastra-Liedahl emission spectrum for optically-thin 
plasmas supplied with XSPEC. The model includes 
both thermal bremsstrahlung and line emission. We 
have assumed the clusters to have 0.3 solar abun- 
dance (Mushotzky et al. 1996). We then weight the 
electron temperature within a given zone (which, of 



course, is given directly by the hydrodynamical evolu- 
tion) by the volume emissivity within that same zone. 
The resulting product, EyT e , is then integrated along 
each LOS through the simulation volume for a cho- 
sen projection angle creating an image which is then 
divided by the total emissivity along each LOS. The 
result is a projected, emissivity-weighted temperature 
map of the cluster at approximately the same res- 
olution of the simulation itself (i.e., 50 kpc). The 
isothermal temperatures used in the following anal- 
ysis are determined by averaging the above temper- 
ature map within a given aperture centered on the 
emission peak. In §|5], we will explore the effect of 
varying aperture size. Similarly, the temperature pro- 
files are generated by an azimuthal averaging of the 
temperature map within concentric annuli centered 
on the emission peak. In both cases, it is important 
to note that in the absence of spherical symmetry, 
temperature structure along the LOS is lost. 

3.2.2. Observed Density Profile 

The procedure used to determine the effective elec- 
tron density measured by observers, n e (r), is some- 
what more complicated, and we simulate the proce- 
dure used by most observers. As before, ey is calcu- 
lated at a given epoch. A LOS integration through 
the simulation volume at a chosen projection angle 
results in an X-ray surface brightness image. The X- 
ray surface brightness profile, S x (r) , is generated by 
an azimuthal averaging of the surface brightness im- 
age within annular rings centered on the peak surface 
brightness. The resulting profile is then fit with an 
isothermal /3-model (Sarazin, 1986), 



S x {r) = S 



, r 
1+ I — 

r r 



-3/3+i 



(5) 



where the core radius, r c , and (3, are free parameters. 
We choose 1 Mpc as a reasonable extent for the ob- 
servable X-ray emission. It is roughly three times the 
initial core radius. Generally, one can not expect a 
good fit to the /3-model much beyond 2-3 core radii. 
It should be noted that the fit parameters, r c and /3, 
can depend strongly on the maximum radius used in 
the fit, and that this will contribute to some of the 
error in H Q . However, in general, r c and (3 are not 
fit independently and vary such that the line integral 
of the /3-model remains somewhat constant. This is 
important since the SZ effect depends on the line in- 
tegral not on the individual parameters. 
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Given the fit parameters, r c and j3, and the peak X- 
ray surface brightness, S , it is possible to determine 
the central electron density, n , using the following 
relation based on Eq. 3 in Henry & Henriksen (1986), 



S Q = N{E 1 /kT e ,E 2 /kT e ) 



£(3/3 - 1/2) 



r c nl{kT e f/\ 



(6) 

Here, N is a normalizing factor dependent on the 
bandwidth [E\,E2) of the X-ray observation, T e is 
the electron temperature, k is Boltzman's constant, 
and r is the complete gamma function. Knowing r c , 
(3, and n , we can now calculate the electron density 
profile as follows, 

-3/3/2 



n e (r) = n a 



1 



(7) 



As with the temperature data, in the absence of 
spherical symmetry, all LOS density inhomogeneities 
are lost. Furthermore, the perceived shape of the den- 
sity distribution (r c , (3, and cllipticity) may differ sig- 
nificantly from the true distribution. 

Since we do not include errors associated with pho- 
ton statistics, all uncertainty in the fit parameters and 
ultimately in the electron density result from inhomo- 
geneities and asphericity in both the density and tem- 
perature distributions or from a real deviation from 
a King model. For consistency in the analysis, the 
X-ray profile is always centered on the highest peak 
in the X-ray emission. The "best" fit parameters are 
determined objectively by a \ 2 minimization routine. 

3.3. An Isolated Cluster: A Test of the Model 

The above analysis was applied to a single, hydro- 
static cluster which is similar to the primary cluster in 
each merger. The test cluster was placed on the same 
grid used in all the merger simulations only now it 
was allowed to evolve in isolation for a time period 
comparable to that of the mergers, ^8 Gyrs. This 
test serves several purposes. First, it demonstrates 
the degree of dynamical stability in our model clus- 
ter. The cluster begins isothermal within ^4r c and in 
hydrostatic equilibrium. Left on its own, it should re- 
main so. However, due to resolution effects (largely in 
the N-body distribution), and infalling ambient mat- 
ter, we do expect some minor evolution of the clus- 
ter. Second, we can test the observational analysis 



described above (£3.2.) by attempting to reproduce 
the temperature and gas density profiles of the clus- 
ter. Finally (and most importantly), in light of the 



inevitable, albeit minor, evolution of the isolated clus- 
ter and the limited observational information, we can 
attempt to accurately reproduce H using the SZ ef- 
fect. This final test allows us to assess the degree 
to which the observed H discrepancy is attributable 
to recent merger activity and the evolving dynamical 
state of the clusters. 

The most noticeable evolution occurs in the ICM 
temperature. Initially, the cluster is isothermal (6.5 
keV) within the central 1.2 Mpc or 4 core radii. After 
4 Gyrs, the time at which core passage occurs in the 
merger simulations, the temperature within this same 
region has risen to 6.9 keV with a zone-to-zone dis- 
persion of 0.22 keV or roughly 3% of the mean. After 
8 Gyrs, the temperature has risen to nearly 7.5 keV 
with a dispersion of 0.29 keV or <4% of the mean. 
Since the ICM heating is rather uniform across the 
cluster, we believe that the primary cause is a general, 
although slight, contraction of the cluster mass distri- 
bution. There is some evidence from the fitting of (3 
models to the synthetic X-ray images that the clus- 
ter core has contracted by ^5% in radius, or ~15% in 
volume. In addition, some heating undoubtedly arises 
from random fluctuations in the gravitational poten- 
tial, particularly in the core. This serves to heat the 
gas by inducing collisions through random motions. 
It is this process that likely produces most of the dis- 
persion about the mean temperature. In any case, 
the cluster remains largely isothermal as indicated by 
the low temperature dispersions, and evolution of the 
mean temperature is small compared to that seen in 
the merger simulations. Finally, Figure shows the 
gas density profile at 0, 4 and 8 Gyrs. As with tem- 
perature, the density evolution is minimal. A slight 
increase in the core density is noted. This is likely re- 
lated to the contraction of the cluster discussed above. 
The /3-model fitting indicates that the cluster remains 
spherical to better than a few percent. 

In order to test the significance of the observed 
evolution to the SZ analysis, we hav e pe rformed 
the observational analysis described in j j3.2.| At each 
epoch, we calculate the true SZ-based H at 100 
random LOS through the cluster core while simul- 
taneously calculating the observed value of H . We 
find at Gyrs, < h obs /h tr ue >= 1-00 ± 0.004, at 4 
Gyrs, < h b s /htrue >— 1-00 ± 0.01, and at 8 Gyrs, 
< hobs/htrue > = 1-03 ± 0.02. In each analysis, the 
cluster was assumed to be isothermal. In the last in- 
stance, we find a +3% discrepancy in the observed 
value of H Q . Our analysis shows that the low level 
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of non-isothermality and small degree of aspheric- 
ity discussed above contribute about equally to this 
systematic effect. We performed the same analysis 
using the projected, emission-weighted temperature 
profiles, and found < h b s /ht rue >= 1.00 ± 0.004, 
0.99±0.02, and 1.00±0.02 for 0, 4 and 8 Gyrs, respec- 
tively. This test demonstrates that 1) our analysis 
procedure can accurately reproduce the temperature 
and density profiles, and 2) the minimal evolution ob- 
served in the test cluster does not significantly influ- 
ence the SZ analysis, and therefore, the results de- 
scribed in §||] can be attributed largely to the post- 
merger evolution. 

4. A DESCRIPTION OF A TYPICAL CLUSTER- 
SUBCLUSTER MERGER EVENT 

A detailed description of the physical evolution of 
cluster-subcluster mergers can be found in Roettiger 
et al. (1997). Here, we will only briefly review the evo- 
lution with an emphasis on those aspects most perti- 
nent to the SZ analysis. 

As stated above, we begin with two, essentially 
isothermal, spherical, hydrostatic clusters which merge 
under the influence of their mutual gravity. The clus- 
ters are composed of two components, one collisional 
(ICM) and one non-collisional (dark matter). Al- 
though they do interact via gravity, the evolution 
of these components differ significantly, particularly 
during the early stages of the merger. We begin 
with the dark matter, or particle, evolution since it is 
the most straight-forward. In the higher mass ratio 
mergers such as those discussed here, the subcluster 
is elongated by tidal forces as it is drawn into pri- 
mary cluster. When coincident with the primary core, 
the subcluster is compressed creating an extremely 
concentrated mass distribution and very deep grav- 
itational potential well. Upon exiting the core, the 
subcluster is largely disrupted appearing downstream 
as a spray of particles some of which are no longer 
bound to the system. Those particles that do remain 
bound return on radial orbits through the remnant 
core. The result is a sustained anisotropy in the par- 
ticle velocity dispersion that supports an elongation 
of the dark matter distribution. As noted by both 
Roettiger et al. (1997) and van Haarlem & van de 
Weygaert (1993), the elongation of the dark matter 
distribution, and, consequently, the gravitational po- 
tential, is long-lived such that the alignment of a clus- 
ter merely reflects the axis of the most recent merger. 



Ultimately the ICM will come into equilibrium within 
the elongated potential creating a prolate gas distri- 
bution. As we will see in §^J, this has strong implica- 
tions for the SZ analysis. 

The evolution of the gaseous components is some- 
what more complicated. As the subcluster begins to 
impinge on the primary cluster, its leading edge is 
compressed and swept upstream as subcluster gas is 
stripped away at large radii. Along the leading edge, 
a bow shock is formed as the two cores begin to merge 
which further compresses and heats the gas. At the 
time of core passage, the gas distribution is strongly 
peaked, as is the gravitational potential at this time, 
and slightly elongated perpendicular to the merger 
axis (Figure |^a from Merger 3; Table 1). As the sub- 
cluster exits the core, the initial bow shock proceeds 
downstream and can be seen as an arc of hot gas. Be- 
hind the shock is a region of rarefaction and adiabatic 
expansion containing much cooler gas (left edge of 
Figure ||d). As noted above, the gravitational poten- 
tial reaches a peak minimum during core passage only 
to rapidly return to near premerger values as the sub- 
cluster exits the core. While near its minimum, gas is 
drawn in from the outer regions of the cluster only to 
be rapidly expelled from the core. The expelled gas 
interacts with residual infall from the subcluster cre- 
ating a second shock that propagates upstream (far 
rig ht Figure |d). 

The fate of the subcluster gas depends on the rel- 
ative cluster parameters. In general, most of the sub- 
cluster gas is stripped before core passage, remain- 
ing as relatively cool gas on the upstream side of the 
merger remnant. Gas, along the leading edge of the 
subcluster is severely heated (Figure ||b) before be- 
ing stripped whereupon it is transported along the 
bow shock and deposited in the outer parts of the 
cluster. In those instances where the subcluster and 
primary cluster have comparable ram pressure, some 
subcluster gas will penetrate and be carried beyond 
the primary core. Since the subcluster gas is no longer 
bound within its own gravitational potential nor is it 
pressure confined by the primary ICM at large radii, 
it eventually expands and cools before falling back 
into the primary where it is subsequently re-heated 
and slowly mixes via random gas motions in the rem- 
nant (Figure ||f,h). All of these effects, the shocks, 
adiabatic expansion, remnant subcluster gas, etc, re- 
sult in a very inhomogeneous, (i.e., non- isothermal) 
temperature morphology which, as we will see in §pj 
has implications for the SZ analysis. 
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5. RESULTS 

The origin of the possible discrepancy in H , at 
least the component associated with the X-ray ob- 
servations, has two primary sources, those being the 
uncertainty in the true temperature and density dis- 
tributions. Here, we attempt to separate and quan- 
tify these two effects by first presenting an analysis 
using the "observed" temperature data, both a single 
isothermal temperature and a high-resolution temper- 
ature profile combined with the true LOS density dis- 
tribution. We then present the same analysis using 
the observed temperature data along with an obser- 
vationally derived density profile. Since the SZ effect 
is highly dependent on the particular LOS through 
the cluster, we present our results as a function of 
projection angle with respect to the observer. And, 
since the systems are evolving, we also examine two 
epochs during the post-merger evolution. 

5.1. Temperature 

Figures |]a-f show the discrepancy between the ob- 
served and true value of H that is due entirely to 
uncertainty in the LOS temperature distribution as 
a function of projection angle. The projection angle 
is defined relative to the merger axis such that +90°, 
0°, and -90° correspond to the observer looking up the 
merger axis (subcluster moves toward the observer), 
perpendicular to, and down the merger axis (subclus- 
ter moves away from the observer), respectively. We 
chose the merger axis as a reference because it is the 
merger remnant's primary axis of symmetry. These 
figures were generated by calculating the ratio of the 
observed to the true H (Eq. ||) along 100 random 
LOS through the cluster core. This results in a dis- 
tribution of points to which we fit a parabola. For the 
sake of clarity, we only plot the parabolic fit. Typi- 
cally, the maximum residual to the fit is ~ ±0.10. The 
solid line represents the resulting discrepancy when 
using a projected, emission-weighted isothermal tem- 
perature estimate taken within a 0.75 Mpc aperture. 
Realistically, we can now expect more detailed tem- 
perature information on a significant number of rela- 
tively nearby clusters (Markevitch 1996). To address 
this possibility, we have performed the SZ analysis us- 
ing azimuthally averaged temperature profiles binned 
by 50 kpc (dashed line). 

There are several sources for the scatter in the 
above relations. Ultimately, they are all tied to the 
extreme sensitivity of the denominator of Eq. [| to 



the particular LOS. As the merger remnant evolves, 
the extreme symmetry seen in the early merger stages 
breaks down giving way to a more irregular temper- 
ature distribution. Consequently, any small variation 
in the LOS can cause a significant change in the line 
integral. On the other hand, the observed temper- 
atures are relatively unaffected by projection since 
they are averaged over large regions of the cluster. As 
stated above, the typical maximum residual to the fit 
is ±0.10, this corresponds to <10% variation in the 
line integral from one LOS to the next. Figures ||a,b 
show the LOS variation in h b s / 'htrue for one merger. 
In Figure ||a, we have chosen LOS in constant lati- 
tude (± to merger axis) while scanning in longitude 
(|| to merger axis). In Figure ||b, we hold longitude 
constant and scan in latitude. Each point represents 
a LOS. 

In order to check for the affects of aperture size on 
the isothermal temperature, we view mergers 4 and 5 
at one projection (-30°, the isothermal temperature is 
not very strongly dependent on projection angle) and 
calculate h b s / htrue for various aperture sizes ranging 
from 0.25 to 1.25 Mpc. The results of this analysis 
are shown in Figure q. We see a general trend of de- 
creasing h i, s / htrue with increasing aperture size. We 
also note the potential for overestimation of H when 
using relatively small apertures. Inagaki et al. 1995 



observed a similar effect (see §5.6.). In the presense of 
a strong temperature gradient, a small aperture will 
result in an overestimation of the cluster temperature 
which can lead to an overestimation of H by 10% or 
more. 

Examination of these figures shows 1) H is sys- 
tematically underestimated in all mergers at all pro- 
jection angles. Using an isothermal temperature es- 
timate, we find errors range from ~ 5-30%. Using 
a temperature profile, errors range from ^5-20%. 2) 
The systematic underestimation of H tends to be 
greater in mergers where the subcluster contains more 
gas relative to the primary, Figures ||a-c,f. 3) Esti- 
mates of H , in all but one case, improve when the 
temperature profile is used, generally by ~10%. The 
improvement is most noticeable in the more disrup- 
tive mergers, (i.e., those where the ratio of the cen- 
tral gas densities is low). Little, if any improvement 
is noted in the very high gas density ratio mergers, 
but these showed relatively small systematic errors to 
begin with. 4) Generally, the error is relatively flat as 
a function of projection angle (latitude), varying by 
<0.1 in h b a /h 

true • 
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5.2. Density and Temperature 

Figure 0a-f shows the discrepancy in H Q when 
both the density and temperatures are observation- 
ally determined. As before, the solid line indicates an 
isothermal temperature was used, a dashed line indi- 
cates a temperature profile. Comparison with Figure 
|^,-f reveals the relative contribution of each to the 
total discrepancy. Several conclusions can be drawn. 
1) When viewing along the merger axis (±90°), the 
discrepancy increases causing a further underestima- 
tion of H , up to 35% is some cases. 2) When viewing 
perpendicular to the merger axis (0°), the discrepancy 
is the same as or less than that due to temperature 
alone, usually <20%. The combination of (1) and 
(2) create a signature by which the H discrepancy is 
maximum when viewed along the merger axis and at 
a minimum when viewed perpendicular to the merger 
axis. The other trends mentioned above still hold. In 
general, the mergers with lower central density ratios 
have larger systematic errors, and they benefit most 
from using a high-resolution temperature profile. 

5.3. Evolution 

We have chosen to perform the SZ analysis at ~2.5 
and 5.0 Gyrs after core passage. The results can be 
seen in Figures ||a,b. There were several reasons for 
selecting these particular epochs. First, as noted in 
Rocttiger et al. (1997), much of the internal gas dy- 
namics have diminished by 2.5 Gyrs, particularly bulk 
flows. The core crossing time for these clusters is < 1 
Gyr which is the absolute lower limit for post-merger 
relaxation of the cluster. The actual relaxation time is 
considerably longer owing, at least in part, to contin- 
ual infalling of subcluster particles and gas for nearly 
2 Gyrs. By waiting until 2.5 Gyrs, we can neglect the 
kinematic SZ effect (§5.4.) in this portion of the anal- 
ysis. A second motivation for the selection of these 
particular epochs is the ill-constrained cluster merger 
rate. Both cosmological numerical simulations and 
observational evidence supports a hierarchical uni- 
verse in which largescale structures form and grow via 
mergers with smaller systems. An interesting ques- 
tion then pertains to the rate at which these mergers 
occur. Unfortunately, ascertaining the true merger 
rate is very difficult. Based on cooling flow evolution 
analysis, Edge et al. (1992) estimate mergers occur 
every 2-4 Gyrs. Observational surveys (both optical 
and X-ray) show substructure in 30-70% of clusters 
(Bird 1993; Mohr et al. 1993; Davis 1994) indicating 



that a large fraction of clusters have undergone signif- 
icant dynamical evolution in the last 2-3 Gyrs. A the- 
oretical analysis by Richstone, Loeb, & Turner (1992) 
shows that in an H, —l universe, only 40% of present 
day clusters are formed at 80% of the Hubble time, in- 
dicating significant accretion has occurred within the 
last 2-4 Gyrs. In an fl o =0.2 universe, about 80% of 
present day clusters are formed by this time. There- 
fore, we chose 2.5 and 5.0 Gyrs to bracket the time 
expected between significant mergers. 

Figure |^a reveals the lack of evolution in Merger 4 
(4:1 mass ratio) between 2.5 and 5 Gyrs. As above, we 
see the signature of a prolate gas distribution, i.e., a 
severe (~30%) underestimation of H at ±90 projec- 
tion. This does not appear to change with time indi- 
cating a sustained anisotropy in the particle velocity 
dispersion which supports the elongated mass distri- 
bution. This would also seem to indicate very little 
evolution in the temperature distribution though only 
a third of the discrepancy was attributable to the tem- 
perature at ±90°. Figure ^3 shows the evolution of 
the H discrepancy in Merger 6 (8:1 mass ratio). At 
2.5 Gyrs, we see the signature of the prolate gas distri- 
bution. However, at 5.0 Gyrs the signature diminishes 
noticeably indicating a more rapid relaxation of the 
gas distribution within a largely spherical potential. 
Once again, there is little evolution in the tempera- 
ture distribution, which contributes the bulk of the 
discrepancy in this case. 

The case for limited evolution is further supported 
by the values in Table 2. Here, we present the mean 
ratio of h (, s /ht r ue averaged over all angles at 2.5 and 
5.0 Gyrs. In each merger, the mean ratio is less 
than unity indicating a systematic underestimation of 
H . However, there is no clear trend regarding the de- 
gree of evolution, except that mergers involving sub- 
clusters with relatively low central gas densities ap- 
pear more likely to show significant positive evolution 
(Mergers 4 and 6 in Table 2). That is, the observed 
discrepancy in H decreases with time. 

5.4. Gas Dynamics 

In addition to the thermal SZ effect, there is also 
a kinematic SZ effect which for small values is an ad- 
ditive term of the form (ctt/c) J n e v pec dl where v pec 
is the velocity of the scattering plasma relative to the 
observer (Sunyaev & Zeldovich 1980). Therefore, the 
estimate of H may be influenced by bulk plasma mo- 
tions induced along the observer's LOS during the 
merger . Roettigcr et al. (1997) showed that merg- 
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ers can generate bulk flows in excess of 1000 km 
s _1 through the cluster core during the early stages 
of a merger. This is particularly true in the cases 
where the subcluster penetrates the primary before 
being stripped of its ICM. Observationally, such flows 
have not been confirmed, although this situation may 
change with the launch of the high resolution X-ray 
spectrometer, Astro-E. 

Unlike the thermal SZ effect, the kinematic effect 
is directional in that it depends on the direction of the 
cluster's motion relative to the observer. Therefore, 
one would not expect to observe a systematic effect 
when studying a large sample of clusters or individ- 
ual clusters where the gas dynamics along the LOS are 
essentially random. However, the kinematic SZ effect 
can lead to an observational bias. Clusters with bulk 
flows directed away from the observer will have sys- 
tematically larger values of AT/Tcmb making them 
easier to detect. This effect is further enhanced by 
the tendency to select hotter, more luminous clusters 
(L x ocT 5 / 2 ). Significant mergers will tend to increase 
the cluster's temperature, and consequently increase 
the X-ray luminosity making them more likely to be 
selected in this sample. Furthermore, clusters merg- 
ing along the LOS will be viewed as end-on prolate 
structures which may be preferentially selected be- 
cause of their enhanced X-ray surface brightness. 

We address the kinematic SZ effect in merging clus- 
ters of galaxies by including an additive term in the 
denominator of Eq. || such that, 
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where v pec is the peculiar velocity of a gas parcel along 
the LOS in the rest frame of the primary cluster's cen- 
ter of mass. Positive v pec is defined to be motion away 
from the observer. We then performed the same anal- 
ysis as above on several selected mergers and epochs. 
In general, we find the kinematic effect is not signifi- 
cant. This is not surprising, since the bulk flow must 
not only be high velocity, but coincident with high 
density. Only in the early stages of mergers 1 and 2 
(those with low ratios of central gas density) do we 
find a noticeable effect, 5-6%, when observing within 
15° of the merger axis and within 1 Gyr of core pas- 
sage. 



5.5. A Merger with a Non-Zero Impact Pa- 
rameter 

We now examine the remnant of a non-zero im- 
pact parameter merger. Although cosmological sim- 
ulations of largescale structure formation would seem 
to indicate that most mergers are largely head-on 
owing to the infall of matter along radial filaments, 
the potential for offaxis mergers does exist. The un- 
usual temperature and density substructure observed 
in A754 (Henry and Briel 1995) is believed to be the 
result of a high angular momentum merger (Henrick- 
sen & Markevitch 1996). The significance of offaxis 
mergers to this study is that residual angular momen- 
tum in the remnant may result in an oblate rather 
than a prolate gas distribution. 

The initial parameters for Merger 7, our offaxis 
merger, can be seen in Table 1. Although a somewhat 
lower mass ratio, Merger 7 is quite similar to Merger 
1. The most significant difference between these two 
mergers is that in Merger 7 the subcluster is given an 
initial velocity of 150 km s _1 perpendicular to the line 
of centers. This is in addition to the 300 km s" 1 along 
the line of centers which was used in all the merger 
simulations. As a result, at the time of closest ap- 
proach the respective centers of mass are separated 
by ~160 kpc or ^0.5r c . This is enough to have a 
significant affect on the early evolution of the merger 
remnant. We find that as the cores interact, there is a 
compression and shear which produces a bar of X-ray 
emission at the interface between the primary and 
subcluster gas components. As the subcluster dark 
matter swings around the primary core, it drags ICM 
with it creating an asymmetric extension in the X-ray 
surface brightness distribution. Near the leading edge 
of the extension, there is an arc of shock heated gas. 
Early in the merger evolution, we find that the dis- 
tinctive X-ray and temperature morphologies and the 
displacement found between the dark matter and gas 
distributions are quite similar to the observed prop- 
erties of A754 (Zabludoff & Zaritsky 1995; Henriksen 
& Markevitch 1996). Unlike the head-on merger, the 
subcluster is not severely disrupted during the initial 
passage. Eventually it reaches turn around and falls 
back into the core of the primary on a nearly radial or- 
bit. The second core passage occurs at approximately 
2.5 Gyrs after the first. During the first passage an- 
gular momentum is transferred to the gas distribution 
while some is lost as dark matter particles are shed. 
Although the angular momentum in the gas should 
dissipate quickly there is still some evidence of rota- 
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tion at 5.0 Gyrs after the initial core passage. There 
are rotational velocities of order 600 km s _1 at radii 
of 1.5-3?' c in the plane of the merger. Velocities are 
lower and less ordered within ~1.5r c . 

Our SZ error analysis is complicated by the asym- 
metric evolution of the off axis merger. As mentioned 
above, the second core passage occurs at approxi- 
mately 2.5 Gyrs after closest approach. Consequently, 
the cluster morphology is still significantly disrupted 
at this time, and it is certainly not recommended that 
such a cluster be included in an SZ-H study. Still, 
for the sake of completeness and consistency, we per- 
form the above analysis at 2.5 Gyrs. We find that the 
results are quite consistent with those of the head- 
on merger. This is as expected since the remnant is 
still largely prolate. There is, however, some oblate- 
ness which causes an overestimation of H by as much 
as 20% when the cluster is viewed within 15-30° of 
face-on (i.e., perpendicular to the oblateness). This 
is caused by an overestimation of the LOS pressure 
integral. Overestimations of both temperature and 
density contribute to this effect. The gas distribu- 
tion that is projected into the plane of the sky ap- 
pears to be more extended than it actually is along 
the observer's LOS, and the temperature is less con- 
taminated by cooler gas along the LOS and therefore 
appears somewhat higher. The analysis was also per- 
formed at 5.0 Gyrs after closest approach. Here, the 
results are even more similar to those of Merger 1 
although the residual oblateness does reduce the un- 
derestimation of H by 5-10% within about 20° of 
face-on. This does not however significantly affect 
the H discrepancy when averaged over all viewing 
angles, compare Merger 7 at 5 Gyrs with Merger 1 at 
5 Gyrs in Table 2. 

5.6. A Comparison with a Previous Study 

Inagaki et al. (1995) also addresses the reliabil- 
ity of H derived from the SZ effect. Our study dif- 
fers from theirs in several significant ways. First, 
our numerical methodology for computing the evo- 
lution of clusters (§|3]) is based on PPM where as 
they employ the numerical simulations of Suginohara 
(1994) which are based on a smoothed particle hydro- 
dynamics (SPH) code. PPM is better able to accu- 
rately record the development and evolution of shocks 
which can dramatically effect the temperature evolu- 
tion within merging clusters of galaxies. Second, we 
have effectively twice the resolution (50 kpc) over a 
much larger region (2 Mpc wide, extending the length 



of our grid, §|3]). Inagaki et al. report 100 kpc reso- 
lution, but this really only applies to the core of the 
cluster. As density decreases, so does their resolution. 
Third, our approach to the analysis is considerably 
different. Inagaki et al. observed two simulated clus- 
ters evolved from cosmological initial conditions but 
without noting their recent history. We have chosen 
to examine a series of clusters with known initial prop- 
erties (based on the observed properties of clusters) 
and very specific merger histories. In this way, we 
can directly relate the observed H discrepancy to the 
cluster's recent evolution. Finally, we have adopted a 
more systematic approach to presenting the analysis 
in order to elucidate projection effects under specific 
observational conditions. The Inagaki et al. study 
was directed toward a more statistical analysis which 
would be applicable to a large sample of SZ-based 
H observations. Our study is meant to show the 
magnitude of errors attributable to any given obser- 
vation. 

With this said, we feel our results are quite con- 
sistent with, although somewhat more extreme than, 
those of Inagaki et al. . They report that non- 
isothermality of clusters, particularly temperature 
gradients, is the most significant source of error in 
H , resulting in as much as a 20% underestimation. 
We find up to 30% in our merger simulations though 
<20% was more typical. They also note that pro- 
lateness of the cluster will cause an underestimation 
of H when viewing the cluster along the major axis. 
However, because of the formalism they choose for 
calculating the core radius, they attribute the under- 
estimation of H to an overestimation of the core ra- 
dius. We, on the other hand, find that fitting a core 
radius and (3 to the X-ray surface brightness of a pro- 
late gas distribution viewed end-on effectively results 
in an underestimation of the LOS extent of the gas 
distribution. Combined with the temperature-based 
underestimation, we find a combined effect of up to 
35%. They quantify the error due to dumpiness and 
asphericity as a <15% overestimation. We also see 
the potential for overestimation of H due to both as- 
phericity and dumpiness. We find viewing a prolate 
distribution, perpendicular to the major axis can re- 
sult in a overestimation of 10-15%. In some instances, 
this is just enough to compensate for the underesti- 
mation of H due to incomplete temperature informa- 
tion. 

6. DISCUSSION 
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6.1. Origin of the H D Discrepancy 

The origin of the H Q discrepancy in these dy- 
namically evolving clusters can be summarized sim- 
ply as a lack of true 3-dimcnsional information and 
the volume-weighted observed quantities. The LOS 
projection effects and limited resolution of the X- 
ray observations tend to mitigate extremes of tem- 
perature. As seen from Eq. [j], the SZ effect is 
most influenced by the hottest, densest gas in the 
cluster. And, although the observed temperature is 
emission- weighted, the emissivity is only weakly tem- 
perature dependent. This, combined with the long 
LOS through the cluster, causes the observed tem- 
peratures to be lower than the effective mean tem- 
perature of the scattering plasma. In the presense of 
even a mild temperature gradient, an isothermal tem- 
perature estimate using a large aperture will include 
a large amount of gas at relatively cool temperatures 
compared to the core, thus weighting the mean tem- 
perature downward. Underestimation of the temper- 
ature, according to Eq. |5[ will lead to an underesti- 
mation of H Q . As the aperture size increases, the ob- 
served temperature decreases which further reduces 
the observed value of H (see Figure ||). Lack of 
spherical symmetry, and the extremely local nature 
of inhomogeneities in the temperature distribution 
cause even the temperature profile to differ signifi- 
cantly from the true LOS temperature experienced by 
the CMB photons. This is best demonstrated in Fig- 
ures ||a and |To| a where we show the observed isother- 
mal cluster temperature within 0.75 Mpc, the az- 
imuthally averaged temperature profile, and the true 
temperature along the observed LOS both perpendic- 
ular and parallel to the merger axis from Merger 4 at 
5 Gyrs. Note that the temperature profile appears 
largely isothermal varying by little more than 1 keV 
across the face of the cluster. In stark contrast, the 
true LOS temperature exhibits several local temper- 
ature peaks. 

Similarly, projection effects play an important role 
in the density analysis. Each of the mergers discussed 
here results in a prolate gas distribution whose major 
axis coincides with the merger axis. Both the degree 
of ellipticity and its evolution are a function of the 
relative total mass and the relative central gas densi- 
ties. As the subcluster increases both in gas content 
and total mass relative to the primary, the degree of 
prolateness in the gas distribution increases. Initially, 
the gas distribution is more elongated than the dark 
matter distribution. As the systems come into equi- 



librium, the gas traces the dark matter fairly well. 
In the case of the 4:1 mass ratio mergers both the 
dark matter and gas distribution remain noticeably 
elongated even after 5 Gyrs. In the 8:1 mass ratio 
mergers, both the gas and dark matter relax into es- 
sentially spherical distributions. 

The signature of the prolate gas distribution is eas- 
ily seen in the results of the SZ analysis. When view- 
ing the elongated cluster along the major axis, projec- 
tion effects cause the gas distribution to appear less 
extended and more centrally concentrated than in re- 
ality. Consequently, both core radius and (3 fit to the 
X-ray profile are too small. This causes the central 
gas density in Eq. || to be slightly overestimated, but 
the primary effect is to underestimate the extent of 
the gas distribution and ultimately H . The effect is 
reversed when viewing along the cluster's minor axis. 
Here, the tendency is to overestimate r c and (3 caus- 
ing a slight underestimation of the central gas density, 
but this effect is more than compensated for by the 
overestimation of the extent of the gas distribution 
along the LOS. This leads to an overestimation of 
H , but it is usually not enough to compensate for 
the underestimate resulting from uncertainties in the 
temperature distribution. Figures and [ic| b show 
a comparison of the true LOS density profile (solid 
line) and the observationally derived profiles parallel 
and perpendicular to the merger axis, respectively. 
This effect varies fairly smoothly between these two 
extremes resulting in the distributions plotted in Fig- 
ures 0a-f. 

6.2. Limitations of this Analysis 

There are several potentially important physical 
processes neglected in this study which may affect 
the determination of H by the SZ effect. For ex- 
ample, we do not include radiative cooling or ther- 
mal conduction, both of which may serve to miti- 
gate some of the effects observed here. (Similarly, 
Inagaki et al. (1995) also neglect these processes). 
This is particularly true of thermal conduction. As 
with radiative cooling, thermal conduction is likely 
to be most important in the dense cores of clusters, 
precisely the region that most influences the X-ray 
emission while also strongly influencing the SZ effect. 
If thermal conduction in clusters is characterized by 
the Spitzer (1962) coefficient of thermal conductivity, 
temperature inhomogeneities will likely be erased on 
relatively short time scales, <1 Gyr (Sarazin 1986). 
However, conduction may be very much limited to 
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the core owing to the strong dependence on density 
and therefore radius. Still, the actual level of ther- 
mal conductivity, which may depend on the detailed 
structure and strength of intracluster magnetic fields 
(Tribble 1989), is uncertain by at least two orders of 
magnitude. For this reason, it is very difficult to ad- 
dress conduction with numerical simulations at this 
time simply because the parameter space is so poorly 
constrained. Radiative cooling in the ICM is better 
understood, but it is difficult to gauge the affect on 
this particular analysis. Cooling may serve to weaken 
radial temperature gradients, and it may also result 
in small scale clumping which poses other problems 
for the SZ analysis. We have chosen our parameter 
space such that cooling is not significant. Cooling 
times range from ~5 Gyrs to greater than a Hubble 
time. 

7. Summary 

We have shown that significant systematic errors 
in the SZ-based value of H can result from a com- 
bination of non-isothermality and asphericity in both 
the temperature and density distributions resulting 
from recent dynamical evolution. Together, these fac- 
tors can cause H to be underestimated by as much 
as 35% with more typical values ranging from 10- 
25%. Although there is a potential for overestima- 
tion of H under specific conditions (in particular, 
oblate clusters resulting from offaxis mergers when 
viewed nearly face-on), the tendency to underesti- 
mate is far more prevalent. This effect will be en- 
hanced by various selection effects which may bias 
SZ samples toward clusters with larger (i.e., more 
detectable) AT/Tcmb- For a given set of X-ray 
properties, the observed H Q estimate will decrease as 
AT/Tcmb increases. Similarly, X-ray based selection 
effects will also tend to bias toward an underestima- 
tion of H . Merger remnants tend to be hotter and 
more luminous which increases their likelihood of be- 
ing selected. Remnants of head-on or nearly head- 
on mergers will be prolate. Prolate clusters viewed 
end-on will be preferentially selected by virtue of an 
enhanced X-ray surface brightness. 

We find that detailed temperature profiles do not 
significantly reduce the error in H , except in the more 
extreme merger examples. It is important to note 
however that no X-ray telescope will supply the true 
3-dimcnsional temperature (or density) information, 
so a reliance on spherical symmetry will remain im- 



portant. In order to minimize the effects of dynam- 
ical evolution, a statistical sample of clusters having 
SZ measurements should be compiled. The mem- 
bers of this sample should appear spherically sym- 
metric in both their X-ray surface brightness, and 
temperature distribution. Elongated clusters with 
twisted isophotes and multiple emission peaks should 
be avoided. Also clusters with large (3 discrepan- 
cies (Pfu ^ Pspec, Sarazin, 1986) as well as those 
with f3 spec much less than or much greater than unity 
should be avoided. Anisotropy in the galaxy velocity 
distribution can be the signature of merger activity 
and when projected into the plane of the sky can re- 
sult in a large (3 discrepancy (e.g., A2256; Roettigcr 
et al. 1995). It will be beneficial for future studies to 
focus on nearby cluster samples which are less subject 
to observational selection effects while having more 
detailed X-ray data. Finally, we suggest that the nu- 
merical modeling of specific systems, such as was done 
with A2256 (Roettiger et al. 1995), will aid in the in- 
terpretation of the SZ base determinations of H . 
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Table 1. Cluster Parameters 



Merger 


Cluster 


M(R<3 Mpc) 


n e (R = 0) 


T e (R=0) 


r b c 


J g 


13 d 

"spec 


n e,p/ n t,S 


ID 


ID a 


10 14 M Q 


10" 3 cm- 3 


keV 


kpc 


R < 1 Mpc 


Ratio 


1 


P 


15.0 


5.0 


9.3 


300 


0.12 


0.7 


1 




S 


3.75 


5.0 


3.7 


190 


0.12 


0.7 




2 


P 


15.0 


5.0 


9.3 


300 


0.12 


0.7 


1.7 




s 


3.75 


3.0 


3.7 


190 


0.09 


0.7 




3 


p 


15.0 


9.7 


6.5 


300 


0.12 


1.0 


2.7 




s 


3.75 


3.7 


3.4 


190 


0.09 


0.75 




4 


p 


15.0 


18.0 


6.5 


300 


0.22 


1.0 


4.9 




s 


3.75 


3.7 


3.4 


190 


0.09 


0.75 




5 


p 


15.0 


9.7 


6.5 


300 


0.13 


1.0 


2.6 




s 


1.87 


3.8 


2.0 


165 


0.15 


0.75 




6 


p 


15.0 


9.7 


6.5 


300 


0.13 


1.0 


11.8 




s 


1.87 


0.8 


2.9 


165 


0.08 


0.5 




7 


p 


15.0 


5.0 


9.3 


300 


0.12 


0.75 


1 




s 


6.0 


5.0 


5.0 


200 


0.08 


0.78 





Table 1: a P=primary S=subcluster. 6 Core radius. c gas fraction by mass within 1 Mpc. d (3 spec = imi p a v 2 
e Ratio of central electron densities Primary/Subcluster. 
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Table 2. < h obs /h tr ue >n 



Merger 


Time" (Gyrs) 


Isothermal T M 


T( r )cd 


1 


2.5 


0.72±0.05 


0.82±0.06 




5.0 


0.64±0.05 


0.77±0.04 


2 


2.5 


0.66±0.04 


0.83±0.05 




5.0 


0.69±0.05 


0.81±0.05 


3 


2.5 


0.81±0.11 


0.90±0.09 




5.0 


0.72±0.09 


0.90±0.09 


4 


2.5 


0.87±0.10 


0.88±0.10 




5.0 


0.88±0.10 


0.92±0.08 


5 


2.5 


0.80±0.07 


0.95±0.07 




5.0 


0.72±0.05 


0.82±0.06 


6 


2.5 


0.87±0.06 


0.80±0.07 




5.0 


0.90±0.06 


0.95±0.05 


7 


5.0 


0.65±0.08 


0.79±0.09 



Table 2: The mean ratio of the observed to true H for each of the first 6 mergers in Table 1, averaged over all 
projection angles at 2.5 and 5.0 Gyrs after core passage. Merger 7 was only analyzed at 5.0 Gyrs. a Time since 
core passage in Gyrs. b Ratio computed using an isothermal temperature estimate. c Ratio computed using a high 
resolution temperature profile. d Uses observed density profile based on an isothermal (5 model. 
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FIGURE CAPTIONS 



Fig. |lj: a) The normalized initial total mass density (solid) and gas density (dashed) for the primary cluster in 
mergers 2-6 (See Table 1). The shape of the total density profile is identical for all simulated clusters, both primary 
and subclusters, only the dimensions and total mass scaling change. The gas density profiles change only slightly. 
This cluster has (3=1. Those with lower f} values have flatter gas density profiles for r > r c . b) The temperature 
profile for the same cluster as in (a). The cluster is isothermal for r < 4r c with a slight temperature gradient at 
larger radii. The other clusters have the same basic temperature profile although scaled in accordance with the 
choice of (3 and the total cluster mass. 

Fig. ^ Evolution of the isolated cluster's gas density profile. At 0.0 Gyrs (solid), at 4.0 Gyrs (dotted), and at 
8 Gyrs (dashed). Core passage for the mergers occurs at ~4.0 Gyrs. 

Fig. ^ Contours of a 2-dimensional slice in gas density (Merger 3, Table 1) taken through the cluster core and 
parallel to the merger axis at a) 0.0 Gyrs, c) 1.25 Gyrs, e) 2.5 Gyrs, and g) 5.0 Gyrs. Times are relative to core 
passage. Contours are uniformly spaced in the logarithm and span 2 orders of magnitude. The subcluster entered 
from the right moving to the left. The corresponding slices in gas temperature can be seen in (b),(d),(f), and (h). 
Contours are linear and most are labeled with the corresponding temperature in keV. All contours are spaced by 
2 keV. Each panel is 4 Mpc on a side. 

Fig. £1 The ratio of the "observed" to "true" H D defined by Eq. | using the emission-weighted isothermal 
temperature (solid line) and the projected, emission-weighted temperature profile (dashed line) as a function of 
projection angle along the merger axis (also cluster's major axis). Lines represent a fit to the distribution of 
randomly drawn LOS through the cluster core. Maximum residuals to the fit are ~0.1. In both cases, the exact 
n e distribution is used. A projection angle of 0° implies the observer is looking perpendicular to the merger axis, 
+90° is looking up the merger axis (subcluster moving toward the observer), -90° is looking down the merger axis 
(subcluster moving away from the observer). Each panel a-f corresponds to a merger, 1-6, in Table 1. 

Fig. ^: The ratio of the "observed" to "true" value of H defined by Eq. || using observationally derived 
temperature and density information as a function of projection angle (a) along the merger axis (i.e., latitude) 
at fixed longitude (0°) and (b) perpendicular to the merger axis (i.e., longitude) at fixed latitude (0°) (Merger 
5, Table 1). Each point represents a single measurement. The + designates use of an isothermal temperature. 
The o designates use of an azimuthally averaged temperature profile. The projection angle in (a) is defined as it 
is in Figure The zero point for the projection angle in (b) is arbitrary owing to the relative symmetry about 
the merger axis. Variation of h b s /h t rue from one pointing to the next is a result of local inhomogeneities in the 
temperature and density distributions. 

Fig. |^: The "observed" discrepancy in H resulting from limited temperature information only as a function 
of aperture size for mergers, 4 (+) and 5 (o) at 2.5 Gyrs (solid line) and 5 Gyrs (dashed line). As the aperture 
increases, temperature gradients generated during the merger, cause the observed temperature to decrease. As the 
temperature decreases so does the observed value of H D . 

Fig. fjj: The ratio of the "observed" to "true" value of H defined by Eq. || using observationally derived 
temperature and density information as a function of projection angle along the merger axis (also cluster's major 
axis). The solid line represents the use of an emission- weighted isothermal temperature. The dashed line represents 
the use of a projected, emission-weighted temperature profile. As in Fig. ||, the lines represent fits to the distribution 
of randomly drawn LOS through the cluster core. Maximum residuals to the fit are ~0.1. The projection angle is 
defined as in Fig. || 

Fig. ||: Evolution of the "observed" discrepancy in H D . a) Evolution in Merger 4. Virtually no evolution is 
noted between 2.5 (solid) and 5.0 (dashed) Gyrs. b) Evolution in merger 5. Here, we see the observed discrepancy 
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increase with time. As the merger remnant evolves, it becomes more spherical and so the overestimate of H at 0° 
is reduced, while the underestimation of H caused by the temperature estimate remains essentially the same. 

Fig. |9|: a) A comparison of the isothermal temperature (+) (R<0.75 Mpc), the observed temperature profile 
(o), and the true temperature in a narrow beam along the merger axis (solid line) for Merger 4 after 5 Gyrs. Note 
the observed temperatures are smoother and tend to underestimate the true distribution, b) A comparison of the 
observed density distribution (+) and the true density along the merger axis (solid line). While the peak density 
is slightly overestimated, the width of the observed distribution is severely underestimated. 



Fig. 10: a) A comparison of the isothermal temperature (+) (R<0.75 Mpc), the observed temperature profile 
(o), and the true temperature (solid line) in a narrow beam perpendicular to the merger axis (solid line) for 
Merger 4 after 5 Gyrs. Similar to Figure |j, the observed temperatures are smoother and tend to underestimate 
the true distribution, b) A comparison of the observed density distribution (+) and the true density (solid line) 
perpendicular to the merger axis (solid line). Contrary to Figure |9|, the peak density is slightly underestimated, 
while the width of the distribution is significantly overestimated. 
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Fig. 1.- 
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Fig. 3.- 



24 





0.50^ 
-100 



1.20 r 

1.10- 

1.00 r 
g 0.90 r 

^°0.80- 

0.70 r 

0.60^ 



-50 50 

Projection Angle (degrees) 

(f) 



100 



0.50 E 



0.50 F 
-100 

1.20 

1.10 - 

1.00 r 

§ 0.90 ^ 
^°0.80- 
0.70 r 

0.60 r 

0.50 



-50 50 

Projection Angle (degrees) 

(e) 



100 



-100 



-50 50 

Projection Angle (degrees) 



100 



-100 



-50 50 

Projection Angle (degrees) 



100 



Fig. 4.— 



25 



0.80 
0.75 

0.70 

ID 

< 0.65 
0.60 
0.55 



0.50 





1 










000000°°° 
















+ 


+++ ++*:++ + + + ++ - 


- + 




+++ + 


o 

o 00 o 




o 


++ + +++ + 


oo _ 











-100 



-50 50 

Projection Angle (degrees) 



100 



0.80 
0.75 

0.70 

< 0.65 

0.60 

0.55 
0.50 



o 



a oO<> o 

°<>~ -<> ^o ooo<> 



o vw v o 



oo 



oo 



oo 



o 



50 100 150 

Projection Angle (degrees) 



200 



2G 




27 






Fig. 7.- 



28 




29 



-1000 



-500 500 

Radius (kpc) 



1000 




0.01 Lj , , , , i , , , , i , , , , i , , , , lJ 

-1000 -500 500 1000 

Radius (kpc) 

Fig. 9.— 



30 



10 




